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I.  INTRODUCTION 


An  All  Weather  Short  Range  Air  Defense  System  (A/W  SHORADS)  concept,  which 
utilizes  data  from  an  on-board  strapdown  Inertial  Measurement  Unit  (IMU)  and 
target  state  updates  from  the  ground  to  provide  a  midcourse  guidance  phase,  is 
currently  being  studied.  The  missile  rates  and  accelerations  are  measured  by 
two  low-coat  multisensors  contained  in  the  IMU.,  In  order  to  reduce  the  naviga¬ 
tion  errors  which  accumulate  from  multisensor  error  sources  during  flight,  a 
pre-launch  multisensor  calibration  is  desirable. 

Rockwell  International  [1]  performed  a  study  to  develop  a  calibration  scheme 
for  use  with  multisensors  on  a  fiber  optics  guidance  missile  (FOG-M)  concept. 

This  scheme,  however,  was  developed  for  a  vertical  missile  orientation  at 
calibration  and  was  not  directly  applicable  to  the  A/W  SHORADS  case.  The  recom¬ 
mended  calibration  procedure  from  Reference  1  was,  therefore,  modified  to  perform 
for  a  level,  instead  of  a  vertical,  missile  calibration  orientation  and  to  be  in 
accordance  with  the  A/W  SHORADS  axis  definitions. '  The  resulting  proposed 
calibration  scheme  is  described  in  this  report.  It  relies  heavily  on  Reference  1 
for  nomenclature  and  for  the  framework  within  which  this  procedure  was  developed. 

Section  II  presents  the  multisensor  configuration  for  navigation  and  the 
calibration  reference  axes,  along  with  the  parameters  to  be  calibrated. 

Section  III  presents  the  accelerometer  calibration  measurement  equations 
and  the  iterative  procedure  used  to  calibrate  the  accelerometer  parameters. 
Section  IV  presents  the  gyro  calibration  equations. 

The  proposed  calibration  procedure  was  programmed  into  an  alignment 
subroutine  in  the  six  degree-of-freedom  (6-DOF)  A/W  SHORADS  digital  simulation 
in  order  to  examine  its  performance.  The  Appendix  contains  a  description  of  the 
sequencing  of  operations  for  the  calibration  along  with  a  listing  of  the 
calibration  section  and  data  outputs  for  several  check  cases. 

II.  SYSTEM  CONFIGURATION 

For  the  A/W  SHORADS  case,  as_in  Reference  1,  each  multisensor  is  mounted 
with  a  calibration  rotation  axis  nominally  perpendicular  to  the  multisensor 
spin  axis  S^.  For  the  A/W  SHORADS  case,  however,  the  aniltisensors  are  mounted 
into  the  missile  in  the  orientation  shown  in  Figure  1.  This  is  the  assumed 
navigation  orientation,  with  ■  -90°  and  9r2  a  0°  *nd  provides  a  redundant 
pitch  axis  instead  of  a  redundant  yaw  axis  as  is  the  case  for  FOG-M. 

In  order  to  perform. the  calibration,  measurements  are  made  by  each  multl- 
sensor  at  three  positions:  *  0s,  -90* ,  -180°.  In  addition,  measurements 

are  made  by  each  multisensor  as  it  is  rotated  between  the  0*  and  -180° 
positions.  The  measurements  made  at  the  stationary  orientations  are  used  in 
calibrating  the  accelerometers  and  for  some  of  the  gyro  calibrations.  The 
data  measured  during  the  180*  rotation  is  used  in  the  gyro  calibration  process. 
The  parameters  which  are  to  be  calibrated  are  shown  in  Table  1  (see  Reference  1) . 


TABLE  1.  PARAMETERS  TO  BE  CALIBRATED 


PARAMETER 

DEFINITION 

KA1*  KA2 

ACCELEROMETER  SCALE  FACTOR 

ACCELEROMETER  MISALIGNMENT  ABOUT 

A1  Aa 

THE  SPIN  AXIS 

BA1’  “bI*  BA2*  BB2 

ACCELEROMETER  BIASES 

®l»  \ 

ANGLE  BETWEEN  R  AND  ?  IS  (tt/2-E) 

RADIANS 

KG1*  KG2 

GYRO  SCALE  FACTOR 

6G1’  5G2 

GYRO  MISALIGNMENT  ABOUT  SPIN  AXIS 

°A1*  DB1*  °A2*  DB2 

GYRO  DRIFT  BIASES 

daigsa 

ON-AXIS  G-SENSITIVE  DRIFT 

daigsb 

CROSS-AXIS  G-SENSITIVE  DRIFT 

III.  ACCELEROMETER  CALIBRATION 


A.  Introduction 

This  section  develops  the  equations  used  in  the  accelerometer 
calibration  process  for  the  A/W  SHORADS  concept.  The  orientation  of  the 
multisensor  axes  at  each  of  the  stationary  measurement  positions  will  be  shown 
relative  to  a  reference  axis  set  defined  as  follows.  Let  multisensor  2  be 
at  its  navigation  orientation,  0p2  ■  0°.  The  reference  roll,  pitch  and  yaw 
axes  for  multisensor  2  are  then  given  as 


roll2  -  -S2 

(1) 

PITCH2  *  B2  *  S2  x  A2 

(2) 

yaw2  -  A2  -  r2  x  S 2  , 

(3) 

and  are  shown  in  Figure  2.  Let  multisensor  1  be  at  its  navigation  orientation, 
*  -90°.  The  reference  roll,  pitch  and  yaw  axes  for  multisensor  1  are  then 
given  as 


ROLL1  -  BX  « 

S1  x  A1 

(A) 

PITCH1  -  Ax  « 

RX  x  Sx 

(5) 

YAWi  -  -S-l 

» 

(6) 

and  are  shown  in  Figure 

3. 

Figure  2.  Reference  axis  set  and  multisensor  2 
axes  at  0g2  ■  0®. 
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Figure  3.  Reference  axis  set  and  multlaensor  1 
axes  at  0R^  *  -90°. 


For  multlsensor  2,  the  6^  *  0°  orientation  Is  taken  to  be  a  perfect 
orientation.  Due  to  mechanical  imperfections,  the  multlsensor  is  not  rotated 
by  exactly  90°  from  one  calibration  position  to  another.  This  effect  is 
modeled  by  Ojg  "  “90°  +  024®*”^°°  an<*  ®R2  “  “180°  +  822  ^  -180® .  Another 
mechanical  imperfection  which  ia  considered  is  the  non-orthogonality  of  the 
calibration  rotation  and  the  spin  axes.  The  actual  angle  between  the  two 
axes  is  taken  to  be  (ir/2  -  E2)  radians.  For  multlsensor  1,  the  0R1  ■  -90° 
orientation  is  taken  to  be  the  reference  point  and  rotation  to  0R2  ■  &11*  °° 
and  to  6^  -  -180°  +  8.12  *  “T®00  define  the  error  angles  Bn  and  822 •  The 
angle  between  and  "5 2  is  (ir/2  -  Ej)  radians.  These  misalignment  angles  thus 
define  the  orientation  of  the  actual  sensor  axes  with  respect  to  the  reference 
axis  set. 


B.  Multlsensor  2  Accelerometer  Calibration 

Multlsensor  2  is  a  pitch/yaw  sensor.  Its  orientation  during  navigation 
is  with  0£2  “  0°  and  an  ideal  (reference)  axis  set  for  use  in  calibration  is 
defined  at  this  orientation  as  shown  in  Figure  2.  Since  the  accelerometer  axes 
are  assumed  to  be  misaligned  about  the  spin  axis  by  an  angle  6^2  >  the  actual 
multlsensor  axes  are  located  with  respect  to  the  reference  axes  as  shown  in 
Figure  A,  where  £42  is  positive  for  a  positive  rotation  about  ?2* 


B2 


Figure  4.  Accelerometer  misalignment  with  respect 
to  spin  axis,  0R2  ”  0*. 
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If  one  writes  the  equations  for  measurements  along  ~&i  and  32,  and 
applies  the  swell  angle  approximation,  the  results  are: 

A21  ■  Kg2  0084^2  ^  6jSnKg2  +  *w  Ol  +  6g2  ®2  +  ®A2  (7) 

B21  *  1^2  ®2  ®®4A2  “  CiSHHa,  +  Bg2  *  1^2  ®2  "  4 42  ®1  +  ®B2  •  (8) 

When  aultlsensor  2  Is  rotated  to  632  “  -90* t  the  axes  are  oriented 
as  shown  in  Figure  5.  The  final  orientation  of  the  actual  axes  with  respect 
to  the  reference  axes  can  be  obtained  with  a  sequence  of  four  rotations.  These 
are  as  follows: 


1.  Perform  a  positive  rotation  of  E2  about  the  YAW.  axis  to  form  an 

Intermediate  axis  set  *2,  *2*  z2*  2 

2.  Perform  a  positive  rotation  of  E2  about  the  X2  axis  to  form  a 
second  Intermediate  axis  set  X3,  Y3,  Z3. 

3.  Perform  a  positive  rotation  of  634  about  the  Y3  axis  to  form 
a  third  Intermediate  axis  set  X4,  Y4,  Z4. 

4.  Perform  a  negative  rotation  by  6a2  about  the  Z4  axis  (this 
corresponds  to  a  positive  rotation  about  the  nultlsensor  spin 
axis)  to  form  the  final  orientation  of  the  nultlsensor  axes  for 
®R2  “  “90°. 


Figure  5.  Msltloensor  2  axes  at  0^  *  -90*. 

If  one  next  forms  the  coordinate  transformation  matrices  for  each  of 
these  rotations,  using  the  small  angle  approximation,  the  measurement  equations 
for  this  orientation  can  written  as 

A24  "  “*42  ®3  “  ®2®2  +  ®A2®2  +  ®24®1  +  ®A2 

*24  *  *A2C2  +  6A2G3  "  *2G3  +  ®2G1  +  *B2  •  (10) 
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When  multisensor  2  Is  rotated  to  0R2  •  -180°,  the  axes  are  oriented 
aa  shown  in  Figure  6.  The  final  orientation  of  the  actual  axes  with  respect 
to  the  reference  axis  set  can  be  obtained  in  this  case  with  a  sequence  of 
thfcee  rotations.  These  are  as  follows: 

1.  Perform  a  positive  rotation  of  2E2  about  the  YAW,  axis  to  form  an 

intermediate  axis  set  X2,  Y2,  Z2.  2 

2.  Perform  a  positive  rotation  of  6^  about  the  X,  axis  to  form  a 

second  intermediate  axis  set  X3,  Y3,  Z3.  1 

3.  Perform  a  positive  rotation  of  about  the  X3  axis  (corresponding 
to  a  positive  rotation  of  6^  about  ’S'2)  to  form  the  final 
orientation  of  the  multisensor  axes  for  0R^  «  -180°. 


YAW 

Figure  6.  Multisensor  2  axes  at  9r2  ■  -J00°. 


If  one  next  forms  the  coordinate  transformation  matrices  for  each  of 
these  rotations,  using  the  small  angle  approximation,  the  measurement 
equations  for  the  0R2  -  -180°  orientation  can  be  written  as 


a22  ”  ~kA2g1  +  $A2G2  ~  P22G3  +  BA2  <n) 

B22  -  XA2G2  "  2E2G3  +  6A2g1  +  BB2.  <12) 

The  calibration  equations  for  multisensor  2  are  then  developed  by 
te  measurements  at  0g2  *  0°  and  —180*  and  the  results  are: 

0.5  (A2j  -  A22  “  @22G3^  "  EA2G1  (13) 

0.5  (B22  -  B21)  +  E2G3  -  6a2g1  <14> 

0.5  (A21  +  A22  +  B22g3>  “  «A2G2  -  H 2  <15> 

0.5  (B21  +  B22)  +  E2G3  -  «A2g2  “  ®B2  <16> 


C.  Multisensor  1  Accelerometer  Calibration 

Multisensor  1  is  a  roll/pltch  sensor  with  the  pitch  axis  a  redundant 
axis.  The  sensor  orientation  during  navigation  is  with  6Rl  -  -90°  as  shown 
in  Figure  1  and  a  reference  axis  set  for  use  in  calibration  is  as  shown  in 
Figure  3.  Since  the  two  multisensors  are  calibrated  with  respect  to  different 
reference  axis  sets,  a  transformation  matrix  relating  the  two  reference  sets  is 
necessary.  This  transformation  can  be  developed  from  a  sequence  of  three 
small  angle  rotations  as  follows: 

1.  Perform  a  negative  rotation  of  Yy  about  the  YAWi  axis  to  form  an 
Intermediate  axis  set  X2,  Yj,  I2' 

2.  Perform  a  negative  rotation  of  ¥p  about  the  Y2  axis  to  form  a 
second  intermediate  axis  set  X3,  Y3,  Z3. 

3.  Perform  a  negative  rotation  of  about  the  X-*  axis  to  form  the 
final  configuration,  which  is  the  ROLL2,  PITCH2,  and  YAW2  axis 
set.  The  angles  Yy,  Yp,  4*R  represent  small  misalignment  angles 
between  the  two  reference  axis  sets  and  the  transformation 
developed  from  the  above  sequence  is  represented  by: 

(ROLL  2  \  [  1  -Yy  'Pp  "I  /\  R^LLj_  \ 

PITCH2  I  -  Yy  1  -Yr  I  fXmCHj.  ]  (17) 

YA»2  /  [-PP  Yr  1  J  \  YAW1  / 

This  transformation  is  used  to  relate  measurements  made  in  one  frame  to  the 
other  frame  during  the  calibration  process. 

Also,  in  accordance  with  the  notation  used  in  Reference  1,  the  gravita¬ 
tional  components  along  the  multisensor  1  reference  axes  are  denoted  by  G^,  ,  G 

The  accelerometer  axes  are  assumed  to  be  misaligned  by  <SA1  about  the 
spin  axis  TTi,  so  the  actual  multisensor  axes  for  0r  -  -90°  are  located  with 
respect  to  the  reference  exes  as  shown  in  Figure  7.  The  measurement  equations 
for  and  for  this  orientation,  can  be  written  as: 

(18) 

(19) 


Al  PITCH 


Figure  7.  Accelerometer  misalignment  about  the 
spin  axis  ®Rl  “  "90#* 


A14  “  KA1  G2  -  5A1  g3  +  ®A1 
b14  “  kai  G3  +  6A1  G2  +  BB1  . 


When  multisensor  1  Is  rotated  to  8^3  -  0°,  the  axes  are  oriented  as 
shown  In  Figure  8.  The  final  orientation  of  the  actual  axes  with  respect  to 
the  reference  axes  can  be  obtained  with  a  sequence  of  four  small  angle 
rotations  as  follows: 

1.  Perform  a  positive  rotation  of  E3,  about  the  PiTCHi  axis  to  form 
an  intermediate  axis  set  X2,  ¥2*  z2* 

2.  .Perform  a  negative  rotation  of  E3  about  the  Z2  axis  to  form  a 

second  intermediate  axis  set  X3,  Y3,  Z3. 

3.  Perform  a  positive  rotation  of  633  about  the  X3  axis  to  form  a 

third  intermediate  axis  set  X4,  ,  Z4. 

4.  Perform  a  positive  rotation  of  6^3  about  the  Y4  axis  to  form  a 
final  orientation  of  the  multisensor  axes  for  8^3  “  0° . 


Figure  8.  Multisensor  1  axes  at  8^3  -  0°. 

If  one  now  forms  the  coordinate  transformation  matrix  representing  this 
rotation  sequence,  the  measurement  equations  for  9ri  “  0°  may  be  expressed  as: 

A11  "  KA1  Gj  +  «A1  G3  +  E3  G3  -  hi  g\  +  Bm  (20) 

®11  *  KA1  G3  '  E1  G2  -  B1  gi  -  «ai  Gj  +  BB1  (21) 

When  multisensor  1  is  rotated  to  6^3  “  -180”,  the  axes  are  oriented 
as  shown  in  Figure  9.  The  final  orientation  of  the  actual  axes  with  respect 
to  the  reference  axis  set  can  be  obtained  in  this  case  with  a  sequence  of  four 
small  angle  rotations  as  follows: 

1.  Perform  a  positive  rotation  of  £3  about  the  PITCH 3  axis  to  form 
an  intermediate  axis  set  X£»  Y2»  ^2* 
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2.  Perform  a  positive  rotation  of  E*  about  the  Z2  axis  to  form  a 
second  intermediate  axis  set  X3,  Y3,  Zj . 

3.  Perform  a  positive  rotation  of  Bi2  about  the  X3  axis  to  form  a 

third  intermediate  axis  set  X4,  ,  Z^. 

4.  Perform  a  negative  rotation  of  «A1  about  the  Y4  axis  to  obtain 
the  final  configuration  for  Opi  -  >180°. 


A1 


Figure  9.  Multisensor  1  axes  at  6^  ”  -180°. 

If  one  now  forms  the  coordinate  transformation  matrix  representing  this 
rotation  sequence,  the  measurement  equations  for  8gx  ■  -180°  may  be  expressed 
as: 

A12  -  -KA1  g\  +  «A1  g]  -  EX  gJ  +  812  gJ  +  Ba!  (22) 

b12  “  kA1  +  E1  G2  +  dAl  G1  "  E1  G1  +  ®B1  .  <23> 

The  calibration  equations  for  multisensor  1  are  then  developed  by  using 
the  measurements  at  6^  •  0°  and  -180°  and  the  results  are: 


I 


0.5  (An  -  Al2  +  (8u  +  Bl2)Gj)  -  Kgi  gJ  +  Ex  gJ  *  R^  (gJ  +  Ex  G*)  (24) 
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(25) 


0.5  (B12  -  Bn)  -  G*  -  «xi  Gi 
0.5  (An  +  A12  +  (Bn  -  Si2)gJ)  -  «A1  G3  "  ®A1  (26) 

0.5  <BU  +  B12)  +  Ex  gJ  -  G\  -  Bb1  .  (27) 

D.  Iterative  Calibration  Procedure 

As  in  Reference  1,  the  accelerometer  parameters  will  be  solved  for  by 
use  of  an  iterative  algorithm.  The  data  necessary  for  this  procedure  are  the 
accelerometer  measurements  taken  at  the  stationary  positions  and  the  gravity 
components  G,,  62*  Go  from  a  separate  I MU.  The  iterative  procedure  for 
determining  the  accelerometer  parameters  for  the  A/W  SHORADS  case  is  shown  In 
Table  2.  Ten  Iterations  are  allowed  In  order  that  all  cases  have  sufficient 
time  to  reach  a  steady-state  solution. 

TABLE  2.  A/W  SHORADS  ACCELEROMETER  ITERATIVE  CALIBRATION  PROCEDURE 


IV.  GYRO  CALIBRATION 
A.  Introduction 

This  section  develops  the  equetlons  used  In  the  gyro  calibration 
process  for  the  A/W  SHORADS  concept.  The  gyro  calibration  process  utilises 


measurements  made  at  the  stationary  positions  to  calibrate  the  gyro  drift 
parameters.  The  gyro  scale  factor  Kq  and  the  gyro  misalignment  about  the  spin 
axis,  &G,  are  calibrated  using  integrated  angular  rate  measurements  made  as 
the  multisensors  are  rotated.  The  reference  axes  defined  for  each  multisensor 
in  Section  III.  A.  are  used  in  the  gyro  calibration.  The  error  angles  and 
misalignments  considered  are  as  defined  in  the  accelerometer  sections. 

B.  Multisensor  2  Gyro  Calibration  from  Stationary  Measurements 

The  multisensor  2  gyro  static  measurements  as  presented  in  this 
section  will  correspond  to  the  same  orientation  sequence,  and  the  appropriate 
figures,  presented  in  Section  III.  B.  for  the  accelerometer  static  measurements. 
The  first  orientation  is  for  6^2  “  0  as  shown  in  Figure  2.  '  The  components  of 
earth  rate,  represented  by  fy,  f2  3*  are  measured  along  the  YAW^,  PITCl^, 

and  ROLL ^  reference  axes,  respectively.  The  values  for  ft \ »  Q2»  an<*  JI3  are 
input  from  an  independent  IMU  as  are  the  gravity  components. 

The  measurement  equations  for  the  0^2  “  0°  orientation  are  as  follows: 

WA21  “  "  fil  +  °A2  +  g2  dA2GSB  +  G 1  Da2GSA  (28) 

WB21  *  G2  +  °B2  "  G1  D^gsb  +  G2  DA2GSA  •  (29> 

The  sign  convention  for  the  cross-axis  sensitivity  term,  D^gsb*  can  be 
visualized  by  assuming  a  fictitious  mass  unbalance  along  the  negative  ^  ax^s 
as  shown  in  Figure  10.  From  Figure  10  it^  can  be  seen  that  the  drift  about 
X2  due  to  a  positive  acceleration  along  B2  will  be  given  by 

awA21  “  g2  dA2GSB  »  (3°) 

i.e.,  a  positive  acceleration  along  B2  along  with  a  positive  cross  axis  G- 
sensitlve  drift  term,  D/^gSB’  wil1  reault  in  a  Positive  drift  about  1L.  The 
drift  about  TT2  due  to  a  positive  acceleration  along  7^  is  seen  from  Figure  10 
to  be 

awB21  "  “^l  °A2GSB. 

For  the  multisensor  2  orientation  of  0^2  “  -90°,  the  gyro  measurement 
equations  are 

WA24  "  +  dA2  +  g2  DA2GSB  "  G3  DA2GSA 

WB24  “  G2  +  DB2  +  g3  DA2GSB  "  °2  DA2GSA  .  (33> 


Figure  10.  Multisensor  2  fictitious  mass  unbalance. 


For  the  nultisensor  orientation  of  6^  “  -180°  the  gyro  measurement  equations 
are 


WA22  “  °l  +  dA2  +  G2  dA2GSB  ~  G1  °A2GSA  (34> 

WB22  *  n2  +  °B2  +  G1  dA2GSB  +  G2  DA2GSA  (3S) 

The  G-sensitive  drift  terms  for  multisensor  2  are  calibrated  by  using 
the  measurements  at  *  0®  and  0^  “  -180°.  The  equations  are  as  follows: 

DA2GSA  “  °‘5  <WA21  “  WA22  +  2G1>%  (3« 

DA2GSB  “  G*3  (WB22  “  wB2l)/Gl  (37) 

The  gyro  bias  terms,  and  Dg2  are  calibrated  from  the  measurements  at 
6r2  -  0°.  The  equations  are  as  follows: 

dA2  “  °*5  <WA21  +  WA22>  '  DA2GSBG2  (38> 

®B2  “0.5  (Wg22  +  WB21^  “  G2  “  DA2GSAG2  *  (39) 

C.  Multisensor  1  Gyro  Calibration  from  Stationary  Measurements 


The  multisensor  1  gyro  static  measurements  as  presented  in  this  section 
will  correspond  to  the  same  orientation  sequence,  and  the  appropriate  figures, 
presented  in  Section  III.  C.  for  the  accelerometer  static  measurements.  The  first 
orientation  is  for  0R^  —90°  as  shown  in  Figure  3.  The  components  of  earth 
rate,  represented  by  0^,  Q*,  0  ,  are  measured  along  the  YAWi ,  PITCH^  and  ROLLi 
axes,  respectively.  These  earth  rate  components  are  related  to  those  along  the 
multisensor  2  reference  axes  by  the  transformation  matrix  given  in  equation  17. 


The  measurement  equations  for  the  0R^  *  -90°  orientation  are: 


WA14  *  «2  + 

WB14  “  fi3  + 

For  the  multisensor 
are 


»aii  -  - 

WB11  “  «3 

For  the  multisensor 
are 

WA12  ’  + 

WBi2  “  a\  + 


°A1  +  G3DA1GSB  +  G2DA1GSA  *40* 

°B1  ~  G2DA1GSB  +  G3DA1GSA  .  (4l> 

1  orientation  of  @R^  “0°,  the  gyro  measurement  equations 

+  °A1  +  G3DA1GSB  +  gidA1GSA  (42> 

+  Dbi  -  gJdaigsb  +  G3DA1GSA  •  <43> 

1  orientation  of  0R^  ■  -180®,  the  gyro  measurement  equations 

°A1  +  G3DA1GSB  “  G1DA1GSA  (44) 

°B1  +  GidA1GSB  +  G3°A1GSA  •  (45) 
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The  G-sensitive  drift  terms  for  multisensor  1  are  calibrated  by  using 
the  measurements  at  0^  •  0*  and  -180s.  The  equations  are: 


°A1GSA  "  °*5(Wa11  "  WA12  +  2°l^Gl 
DA1GSB  *  °*5<WB12  "  wB11>/gi  . 


(46) 

(47) 


The  gyro  bias  terms,  D^i  and  are  calibrated  from  the  measurements  at 
0jjl  ■  -90°.  The  equations  are: 


DA1  "  O.SCWau  +  WAi2>  “G3DA1GSB 

®B1  *  g*5(WB11  +  WB12)  -  ^3  -  G3DA1GSA  . 


(48) 

(49) 


0.  Multisensor  2  Gyro  Scale  Factor  and  Misalignment  about  Spin  Axis 
Calibration 

The  calibration  of  Kg2  and  5q2  utilizes  angular  rate  measurements  that 
are  integrated  into  angular  displacement  measurements  as  the  multisensor  is 
rotated  from  0^2  *  0°  to  ®?_2  “  -180°  +  822*  The  angular  rate  measurement 
equations  during  rotation  are: 


WA2  "  KG2GA2  +  6G2  (®R2  +  fi2)  +  DA2  +  dA2GSAA2  +  dA2gsjj®2 

WB2  "  KG2  (®R2  +  "  ^Gl^Al  +  082  “  DA2GSB^2  +  DA2GSAB2 

where  QA2  *  -fljCOSejg  +  J^SIHejg  . 


(50) 

(51) 

(52) 


The  measured  angular  rates  are  then  Integrated  over  the  rotation 
interval  tF2  -  tQ2  ”  t^  to  arrive  at  the  following  angular  displacement 
measurements : 


SF2 

9A2  WA2dt  ”  ^2  0QA2  +  «G2  <A0R2  +  02tR2)  *  DA2CR2  +  DA2GSA  AVA2  +  dA2GSbAVB2 
^02 


(53) 


eB2  “  ^F^B2dt  “  KG2  (A0R2  +  n2cR2>  '  6G2e  fi  A2  +  °B2tR2  +  DA2GSAAVB2  - 
c02 

(54) 


where 

*F2 

0OA2  “  /  nA2dt 

*02 


A0R2  "  /  ^dt  -  (-180°  +  8,2) 
t02  22 


(55) 


(56) 
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(57) 


AVA2  "  /t?2 
*02 


A2<lt 


avb2  “  y* 
*02 


*F2 

B2dt 


(58) 


In  order  to  integrate  A2  and  B2  general  equations  for  the  accelerometer 
measurements  must  be  developed  that  will  satisfy  the  three  stationary 
orientations  and  the  dynamic  orientations  that  occur  during  rotation.  The 
general  equations  for  multisensor  2  are 


a2  "  KA2  <G1  COS(&R2t)  +  G3SIN  (e^t))  +  +  B^  +  E2  ((^SINCe^t)) 

(59) 


B2  "  kA2G2  ■  6A2  (GlCOS(eR2t)  +  G3SIN(6R2t))  +  BB2  -  E2  (G1SIN(6R2t)) 


-E2G3  (1  -COS(0R2t)). 


(60) 


These  correspond  to  measurements  along  the  and  B2  multisensor  measurement 
axes,  respectively. 

After  the  six  second  rotation  and  simultaneous  integration  occur,  all 
the  terms  in  equations  (53)  and  <54)  are  known  from  measurements  or  calculation 
except  for  5q2  and  Kg2*  The  gyro  drift  parameters  were  calibrated  from 
stationary  measurements.  The  error  angle  822  available  from  laboratory 
calibration.  The  earth  rates  are  computed  inputs  from  the  I MU  module.  Thus, 
equations  (53)  and  (54)  can  be  used  to  calibrate  6g2  and  K<;2* 

Matrix  inversion  is  chosen  as  the  most  effective  method  to  solve 
for  6q2  and  Kq2 •  Equations  (53)  and  (54)  are  rewritten  as 

9  A  "  ®A2  DA2tR2  "  dA2GSAAvA2  “  DA2GSBAVB2  "  KG2ef5A2  +  SG2  (A0R2  +  n2tR2) 

(61) 

®B  “  9B2  “  DB2tR2  +  DA2GSBAVA2  ‘  DA2GSAAVB2  "  KG2  <A6R2  +  n2CR2^  “  5C2enA2 

(62) 
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and  can  be  expressed  in  matrix  form  as 


+  fl2 fcR2 


a6R2 


+  n2cR2 

”  6qA2  . 


(63) 


After  the  matrix  inversion  is  performed,  one  obtains 


1 

DETERMINANT 


-0Qa2  “(a0R2  +  ^2fcR2^ 

-(A0r2  +  n2tR2)  0fiA2 


where 

Determinant  -  +  (A0R2  +  £i2tR2)2)  •  (65) 

One  can  now  solve  for  Kq2  and  6q2  as  follows: 


kG2  “  ^0A2  “  °A2tR2  "  dA2GSAaVA2  ~  dA2GSBAvB2>  9nA2  +  (eB2  “  DB2tR2  +  dA2GSBaVA2 
“DA2GSAAVB2^  ^&0R2  +fl2tR2))/((er<A2)2  +  *A0R2  +  ^2tR2)2^ 

6G2  "  ^9A2  “  °A2tR2  ~  DA2GSAAVA2  “  DA2GSBAVB2>  ^0R2  +  fi2tR2)  “  (0B2  ~  DB2tR2 

+dA2GSBAVA2  ~  DA2GSAAVB2>  ^ /<(0«A2> 2  +  (A0R2  +  .  (67> 

B,  Multisensor  1  Gyro  Scale  Factor  and  Misalignment  about  Spin  Axis 

Calibration 

The  calibration  of  Kq2  and  6G2  utilises  angular  rate  measurements  that 
are  integrated  into  angular  displacement  measurements  as  the  multisensor  is 
rotated  from  "  g^  to  0ri  -  -180®  +  8i2.  The  angular  rate  measurement 
equations  during  rotation  are: 

WA1  "  *G1^A1  +  *G1  (0R1  +  n3>  +  °Al  +  °A1GSA  A1  +  DA1GSB  B1  <68> 

WB1  "  *^1  ^®R1  *  _  6G1°A1  +  °B1  _  DA1GSB  A1  +  °A1GSA  B1  ,  *69* 

where  «  -flj  COS  6^  -  °2  SIN  0R1*  (70) 
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The  Measured  angular  rates  are  then  Integrated  over  the  rotation 
Interval  tpi  -  t„.j  “  tjQ  to  arrive  at  the  following  angular  displacement 
measurements: 


tpi 

0A1  "J  WAldt  "  KgI^aj  +  %  <A6R1  +  ^Ri)  +  DAltRl  +  DAlGSAAvAl  +  DA1GSBAVB1 
*01 

(71) 

*Fl 

®B1  »/WBldt  "  *61  (A0R1  +  G3*R1^  '  6GleflAl  +  ^l^Rl  +  DA1GSAAVH  "  DAlGSBAVAl 

(72) 


“01 


where 

*F1 

flrtAi  " 

f  «Aldt 

*01 

-*F1 

&0R1  * 

rr 

1 

CD* 

n 

ft 

(-180°  -  fill  +  8 12 J 


e  tFl 

AVAl  "  J  Aldt 
*01 


avbi  -  /P1Bldt 

toi 


(73) 


(74) 


(75) 


(76) 


in  order  to  Integrate  Ai  and  ,  general  equations  for  the  accelerometer 
measurements  swat  be  developed  that  will  satisfy  conditions  at  the  three 
stationary  orientations  as  well  as  providing  data  during  the  >180*  rotation. 

The  general  equations  for  nultisensor 'l^are 


A1  *  *G1  <G1  008  ($Rlt)  *  G2  SIH  (®gi*))  +  «A!  gJ  (1  +  2  SIN  (dglt)) 

+  ®A1  +  E1  (®J  COS  (^t))  (77) 


19 


mmmmm 


Bi  -  Kai  g\  -  (gJ  COS  (S^t)  +  g\  SIN  (S^t))  +  BB1  -  Ex  (GJ  COS  (0^)) 

-ElGJ  (1  +  SIN  (ftRit))  .  (78) 

These  correspond  to  measurements  along  the  and  Bj  multisensor  measurement 
axes,  respectively. 

After  the  six  second  rotation  and  simultaneous  Integration  occur 
all  the  terms  in  equations  (71)  and  (72)  are  known  from  measurements  or 
calculation  except  for  6q^  and  The  gyro  drift  parameters  were 

calibrated  from  stationary  measurements.  The  error  angles  8n  and  8^2  are 
available  from  laboratory  calibration.  The  earth  rates  are  computed  inputs 
from  the  1MU  module.  Thus ,  equations  (71)  and  (72)  can  be  used  to  calibrate 
«G1  and  K^. 

As  in  Section  IV.  D. .  matrix  inversion  is  used  to  solve  for  6qi  and  K^. 
Equations  (71)  and  (72)  are  rewritten  as: 


®A  “  °A1  “  °A1  CR1  -  DaIGSAAVai  -  Da1GSbAvB1  “  KGle^i  +  «G1  (ABjh  +0  Jty) 


(79) 


0g  “  ®B1  *  Db!  cr1  +  Da!GSBaVa1  -  I>a1GSAAVB1  *  (A0R1  +  ^tR1)  -  «g1^A1  , 


which,  in  matrix  form,  la  expressible  as: 

(A0ri 


/  6a\  9%1 

\«b  j  (A0R1  +  nJtRi) 


+  nJtRi)  / Kgi\ 

-0nAi  \6g1/  . 


After  the  matrix  Inversion  la  performed,  one  obtains: 

(*gi\  x  r-enAi  -(A6ri  +  oatRi)"!  /  0a  \ 

.  J"  DETERMINANT  I  1  I  1  0  1 

®G1/  [-(A8R1  +  OjtRi)  0qa1  J  \  B/  » 


(80) 

(81) 


(82) 
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where 


(83) 


X  2 

DETERMINANT  -  -(®flA1>2  +  <A0R1  +  0^81)  >• 

One  can  thus  solve  for  K<ji  end  Sgl  as  follows: 

*G1  •  ((®A1  "  °A1  tRl  -  DA1GSAAVa1  "  DAlGSBAVBl>  SQaI  +  (eBl  -  ^Bl  fcRl 

~DAlGSAAVBl)  (a®R1  +  ®3fcRl^ ^  f  (  (®5lAl)  ^  +  (A8R1  +  n3fcRl^2^ 

5G1  ■  <(eAl~  °A1  tRl  “  dA1GSAAvA1  -  DAlGSBAvBl)  (A6R1  +  ojtn)  -  (eB1  - 
+DA1GSBAVA1  -  Da!GSAAVB1)  6fiAl)/«eftAl>2  +  <A6Rl  +  0J*B1>2>  • 


+  DA1GSBAVA1 

(84) 

°B1  tRl 

(85) 
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APPENDIX.  CALIBRATION  PROCEDURE 


The  calibration  procedure  described  in  this  report  was  programed  into 
an  alignment  subroutine  in  the  A/W  SHORADS  6 -DOF  digital  siinilation.  The 
nonenciature  used  in  the  subroutine  is  equated  to  the  nomenclature  used  in 
this  report  in  the  program  dictionary  given  in  Table  A-l,  a  listing  of  the 
calibration  algorithm  is  given  in  Table  A-2.  The  results  of  several  check 
runs  are  presented  in  Table  A-3. 

The  procedure  for  calibrating  the  gyro  and  accelerometer  parameters  was 
simulated  in  the  following  sequence: 

(1)  Known  parameters  are  set  (laboratory  calibrations  permit  knowledge 
of  error  angles  and  misalignment  angles  and  literature  from  the  manufacturer  of 
the  multisensor  provides  additional  information  necessary) . 

(2)  The  IMU  provides  measurements  of  earth  and  gravity  rates. 

(3)  Static  measurements  are  made  at  0°,  -90°,  -180°. 

(4)  Angular  rate  measurements  are  made  as  the  multisensor  rotates 

from  0°  to  -180® ,  and  these  measurements  are  integrated  into  angular  displacements 
(the  total  time  necessary  for  the  rotation  is  six  seconds;  the  rate  of  rotation  is 
-30®  per  second) . 

(5)  Ten  iterations  of  the  accelerometer  calibration  equation  set  are 
made  in  order  to  provide  sufficient  time  for  any  extreme  errors  to  settle  out 
to  good  solutions.  From  this  Iterative  procedure  all  the  accelerometer 
parameters  are  calibrated. 

(6)  The  gyro  drift  parameters  are  calibrated  from  static  measurements 
made  by  the  multisensor  and  the  estimated  gravity  terms  from  step  5. 

(7)  All  terms  in  equations  (S3)  and  (54)  for  multisensor  2  and 
equations  (71)  and  (72)  for  multisensor  1  are  known  from  measurements  or 
calculation  except  for  Sgl  end  KgI.  Thus,  these  equations  are  used  to  calibrate 
fiGl  *“d  K^. 

(8)  Calibrated  error  terms  for  both  accelerometer  and  gyro  models  are 
computed. 
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TABLE  A-l.  CALIBRATION  PROGRAM  DICTIONARY 


FORTRAN  NAME 

VARIABLE 

DEFINITION 

A±j  (1-1,2;  j-1,2,4) 

Measured  A|  axis  acceleration 
measurement  at  orientation  j 
(j»l  for  =  0°,  j“2  for  0Ri*“  -180°, 

j-4  for  0Ri  »  -90°) 

ABO 

Go 

Magnitude  of  earth  gravity  at  present 
position. 

AB1,  AB2 ,  AB3 

Gl»  g2,  G3 

Components  of  gravity  (nominally  along 
roll,  pitch,  and  yaw)  corresponding  to 
the  ideal  coordinate  frame  defined  for 
multisensor  2. 

AB1PRM,  AB2PRM 

G1  G1  G1 

Components  of  gravity  corresponding  to 

AB3PRM 

the  ideal  coordinate  frame  defined  for 
multisensor  1. 

ABIASX 

Accelerometer  bias  along  ACS  X-axis 

ABIASY 

Accelerometer  bias  along  ACS  Y-axis 

ABIASZ 

Accelerometer  bias  along  ACS  Z-axis 

AC1 

Pitch-yaw  cross  G-sensitivity 

AC2 

Roll-pitch  cross  G-sensitivity 

AL1 

Pitch-yaw  inline  G-sensitivity 

AL2 

Roll-pitch  inline  G-sensitivity 

ALNFLG 

Align  control  flag 

AMU,  AML2 

fiAl»  6A2 

Accelerometer  misalignments  about  the 
spin  axes. 

AS FI,  ASF2 

KA1»  kA2 

Accelerometer  scale  factor 

ASFX 

Accelerometer  scale  factor  error  along 
ACS  X-axis. 

ASFY 

Accelerometer  scale  factor  error  along 
ACS  Y-axis. 

ASFZ 

Accelerometer  scale  factor  error  along 
ACS  Z-axis. 

Measured  axis  acceleration  measure¬ 

ment  at  orientation  j . 

BAU-BA33 

Elements  of  BCS  to  ACS  transformation 
matrix. 
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TABLE  A-l.  CALIBRATION  PROGRAM  DICTIONARY  (CONTINUED) 


FORTRAN  NAME 

VARIABLE 

DEFINITION 

BA1,  BA2 

bA1.  bA2 

Accelerometer  biases  for  axes  A^. 

BB1,  BB2 

®B1»  bB2 

Accelerometer  biases  for  axes  B^. 

BG11-BG33 

Elements  of  BCS  to  GCS  transformation 
matrix. 

BTA11,  BTA12 

Bll»  e12 

Rotation  angle  errors  for  multisensor 

BTA21,  BTA22 

e21»  S22 

X# 

Rotation  angle  errors  for  multisensor 
2. 

CTHRl,  CTHR2 

a6r1.  a8R2 

Total  mechanical  rotation  angle. 

DAI,  DA2 

Dai.  dA2 

Gyro  drift  bias  for  axes  Ai  and  A2. 

DB1,  DB2 

Db1.  ^2 

Gyro  drift  bias  for  axes  Bi  and  B2. 

DA1GSA,  DA2GSA 

dA1GSA»  dA2GSA 

Gyro  G-sensitlve  drift  sensitivity  to 
acceleration  along  the  angular  rate 
axis. 

DA1GSB,  DA2GSB 

DA1GSB*  DA2GSB 

Gyro  G-sensitive  drift  cross-axis. 

DCMPPT 

YP 

Misalignment  angle  in  pitch. 

DCMPRL 

*R 

Misalignment  angle  in  roll. 

DCMPYA 

*Y 

Misalignment  angle  in  yaw. 

DGSA1,  DGSA2 

Al»  A2 

Multisensor  measurements  along  the  Aj_ 
axis. 

DBSB1,  DGSB2 

®1»  ®2 

Multisensor  measurements  along  the  Bi 
axis. 

DOME1,  DOME2 

dai.  «A2 

Projection  of  earth  rate  onto  the  A^ 

axes. 

El,  E2 

El,  E2 

Non-orthogonality  of  rotation  axis  to 
spin  axis. 

EAA 

e«A 

Misalignment  of  the  pitch-yaw  multi- 
sensor  due  to  temperature  variations. 

Misalignment  of  the  roll-pitch  multi¬ 
sensor  due  to  teapersture  variation a. 


TABLE  A-l.  CALIBRATION  PROGRAM  DICTIONARY  (CONTINUED) 


FORTRAN  NAME 

VARIABLE 

DEFINITION 

EAW1 

Separation  of  ACS  and  GCS  Y-axis  due 
to  temperature  variations. 

EAW2 

Separation  of  ACS  and  GCS  Y-axis  due 
to  temperature  variations. 

EAH3 

Separation  of  ACS  and  GCS  Z-axis  due 
to  temperature  variations. 

EB1 

eBl 

Misalignment  in  azimuth  between  pitch- 
yaw  multlsensor  spin  axis  and  BCS  X- 
axis. 

EB2 

eB2 

Misalignment  in  elevation  between 
pitch-yaw  multisensor  spin  axis  and 

BCS  Y-axis. 

EOA 

eoA 

Orthogonality  error  in  the  pitch-yaw 
multisensor  normal  to  the  spin  axis. 

EOAP 

eoA 

Orthogonality  error  in  the  roll-pitch 
multisensor  normal  to  the  spin  axis. 

ETAXAC 

nxBCS 

Component  of  gravity  along  BCS  X-axis. 

ETAYAC 

nyBCS 

Component  of  gravity  along  BCS  Y-axis. 

ETAZAC 

nzBCS 

Component  of  gravity  along  BCS  Z-axls. 

GBIASP 

Gyro  roll  rate  bias. 

GBIASQ 

Gyro  pitch  rate  bias. 

GBIASR 

Gyro  yaw  rate  bias. 

GL1 

Accelerometer  roll  rate  sensitivity. 

GL2 

Accelerometer  pitch  rate  sensitivity. 

GML1,  GML2 

8G1»  fiG2 

Gyro  misalignments  about  the  spin 
axis. 

GSA1,  GSA2 

AVa1. 

Value  of  Integrated  acceleration  over 
interval  tg^. 

GSB1,  GSB2 

Avbi»  avB2 

Value  of  integrated  acceleration  over 
interval  tgi. 

GSP1,  GSF2 

KG!,  KG2 

Gyro  scale  factor. 
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TABLE  A-l.  CALIBRATION  PROGRAM  DICTIONARY  (CONTINUED) 


FORTRAN  NAME 

VARIABLE 

DEFINITION 

GSFP 

GSFQ 

GSFR 

Gyro  roll  rate  scale  factor  error. 

Gyro  pitch  rate  scale  factor  error. 
Gyro  yaw  rate  scale  factor  error. 

LAST 

Flag  that  controls  six  second 
integration . 

LST 

Flag  that  controls  parameter 
initialisation. 

OME1 ,  OME2 

0^A1»  0fiA2 

Value  of  Integral  of  0^  obtained 
over  interval  tp-f . 

P 

P 

Missile  roll  rate. 

PHI  I 

Misalignment  in  attitude  between 
multisensor  and  BCS  frame. 

PSI12,  PSI13,  PSI21 
PSI23,  PSI31,  PSI32 

Elements  of  BCS  to  GCS  transformation 
matrix  (compensate  for  non-orthog¬ 
onality  and  misalignment  errors) . 

PSII 

*i 

Misalignment  in  azimuth  between 
multisensor  and  BCS  frame. 

Q 

Q 

Missile  pitch  rate. 

QDT 

Quantization  time  rate. 

R 

R 

Missile  yaw  rate. 

RATE 

®Ri 

Mechanical  rotation  rate. 

Rl,  R2,  R3 

tU»  02*  03 

Components  of  earth  rate  (nominally 
along  roll,  pitch,  yaw)  corresponding 
to  the  ideal  coordinate  frame  defined 
for  multisensor  2. 

R1PRM,  R2PRM, 

R3PRM 

aj.  aj.  oj 

Components  of  earth  rate  corresponding 
to  the  ideal  coordinate  frame  defined 
for  multisensor  1. 

RD 

0Ri 

Angular  orientation  of  rotation 

mechanisms . 

Roll*  Pitch*  Yew* 

Ideal  coordinate  frame  defined  for 
multisenaor  2. 
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TABLE  A-l.  CALIBRATION  PROGRAM  DICTIONARY  (CONTINUED) 


FORTRAN  NAME 

VARIABLE 

DEFINITION 

■oil1.  Pitch1, 

Ideal  coordinate  frame  defined  for 

Yawl 

aultlsensor  1. 

SDVX 

Gravity  component  along  the  ACS  X- 
axl8. 

SDVY 

Gravity  component  along  the  ACS  Y- 
axie. 

SDVZ 

Gravity  component  along  the  ACS  Z- 
axis. 

TA1,  TA2 

eAl»  ®A2 

Value  of  Integrated  angular  rates 
over  Intervals  tgj.. 

TB1,  TB2 

®B1»  ®B2 

Value  of  Integrated  angular  rates 
over  intervals  tg*. 

TALIGN 

c01»  c02 

Lover  limit  of  Integration  during 
rotation. 

TF1,  TF2 

‘Pit  *F2 

Upper  limit  of  Integration  during 
rotation. 

THT12 ,THT13 ,THT21 

Elements  of  BCS  to  ACS  transformation 

TBT23 ,THT31 ,THT32 

matrix  (compensate  for  non-orthogo¬ 
nality  and  misalignment  errors) . 

THTAI 

ei 

Misalignment  In  elevation  between 
multisensor  and  BCS  frame. 

TIME 

Simulation  time. 

TR1,  TR2 

tRl*  *R2 

Times  to  rotate  from  0*  to  -180* 
orientation. 

BAjj.(l«l,2jJ*l,2,A) 

Measured  Ai  axis  angular  rates  at 
orientation  j. 

WBlj(l-l*2ij-l,2,4) 

Measured  Bi  axis  angular  rates  at 

orientation  j. 

' 

2i 


'UUOUUU  u  o  o  u 


TABLE  A-2.  CALIBRATION  ALGORITHM 


PROGRAM  LISTING 
IF<LST.  EG  1)00  TO  77 

IF<SNGL(TIME).  LT.  (TALIQN-GDT+O.  0001) )  GO  TO  SOSO 
IFIALNFLO.  GE.  0.  5)00  TO  S04S 
WRITE<6> 99)  8NGL (TIME) 

99  FORMAT <2X/  'ROCKWELL-COLINS  CAL.  CALLED  AT  '.Fll.  4.  ' 
C ****** 

C ******  SET  PARAMETERS  VIA  MCARLO  VALUES 
C***««* 

ASF1  -  ASFY  +  1.0 
ASF2  •  ASFZ  +  1.  0 
BA1  -  ABIASY 
BB1  »  ABIASX 
BA2  -  A8IASZ 
BB2  «  ABIASY 
AML1  »  EAAP 
AML2  *  EAA 
0SF1  -  OSFQ 
G8F2  -  OSFR 
OBI  *  GBIASP 
DAI  -  GB1ASG 
DB2  -  0BIA8G 
DA2  -  0BIA8R 
GML1  -  EAW1 
GML2  -  EAW2 
DA20SA  «  AL1 
0A20SB  -  AC1 
DA108A  -  AL2 
DA 10 SB  -  AC2 


THE  ACTUAL  GRAVITY  AND  EARTH  RATE  COMPONENTS 
ACTINO.  RESPECTIVELY.  ALONG  AND  ABOUT  THE  MULTI SENSOR 
SENSITIVE  AXE8  ARE  AS  FOLLOH8: 


AB3  -  BAI 1 1 *ETAXOO  +  BAI12*ETAY00  +  BAI13*ETAZOO 
AB2  -  BAI21*ETAXOO  +  BAI22*CTAYG0  +  BAI23*ETAZ00 
AB1  -  BAI314ETAX00  +  BAI32*ETAY00  +  BAI33*ETAZ00 
ABO  -  8GRT(AB1*AB1+AB2*AB2+AB3*AB3) 

AB1PRH  -  +DCMPPT*AB3  -  DCMPRL*AB2  +  ABl 

AB2PRM  -  -DCMPYA*AB3  +  AB2  +  DCMPRL*AB 1 

AB3PRM  -  AB3  +  DCMPYA*AB2  -  DCMPPT*AB1 

R3  -  BAI1K0MEGXB  +  BAI  1 2*0MEGYB  +  BAI13*0ME0ZB 

R2  -  BAI21*OMEOXB  +  BAI 22*0ME0YB  +  BAI23*OMEGZB 

R1  -  BAI31*0ME0XB  ♦  BAX32*0ME0YB  +  BAI33*OHEOZB 

R1PRM  >  +DCMPPT*R3  -  DCMPRL*R2  +  R1 

R2PRM  -  -DCMPYA*R3  +  R2  +DCMPRL*R1 

R3PRM  ■  R3  ♦  DCMPYA*R2  -  DCMPPT*R1 

LST  -  LST  ♦  1 

ACCELEROMETER  NOOEL 

OUTPUTS  MEASURED  VALUES  (INPUTS  TO  CALBRA '  EONS) 

A21  -  A8F2*AB 1 +AML2*AB2+BA2 


SEC.  ') 
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A22  -  -ASF2* AB 1 + AML2*AB2-BT A22*AB3+B A2 

A24  -  -ASF2#AB3-E2*AB2+AML2*AB2+BTA24*AB1+BA2 

B21  *  ASF2*AB2-AML2*AB1+BB2 

B22  -  ASF2»AB2-2. 0*E2*AB3+AML2*AB1+BB2 

B24  -  ASF2*AB2+AML2*AB3-E2*AB3+E2*AB1+BB2 

All  -  A8Fl*ABlPRM+AMLl*AB3PRH+El*AB3PRH-BTAll*AB2PRM-*-BAl 

A12  -  -ASF 1 *AB 1PRM+AML 1 *AS3PRH-£ 1 *AB3PRH+BTA1 2*AB2PRM+BA 1 

A14  -  A9F1*AB2PRM-AML1*AB3PRH+BA1 

Bll  -  ASF  1  *AB3PRM-E  1  *AB2PRM— E  1  *AB  1 PRM— AML  1  *AB  1 PRH+-BB 1 
B12  -  A8F1 *AB3PRM+E 1 *AB2PRH+Aftt.l *AB 1PRM— E1*AB 1 PRM+BB 1 
B14  -  ASF 1 *AB3PRM+AML 1 *AB2PRM+BB 1 

C«*M 

C*****  OYRO  MODEL 
C»*** 

UA21  -  -R 1 +DA2+DA2GSB*aB2+DA2CSA*AB 1 
WB21  »  R2+DB2-DA2CSB*AB1+DA20SA*AB2 
UA22  -  Rl+DA2fDA2CSB*AB2-DA20SA*ABl 
MB 22  -  R2+DB2+DA2CSB  *AB  1  ♦DA20SA*AB2 
WA24  *  -R3+DA2+DA20SB*AB2-DA20BA*AB3 
UB24  -  R2+DB2+DA20SB •AB3>DA20SA*AB2 
WA14  ■  R2PRM+DA 1 +DA 1 088 *AB  3PRM+DA  1 08A* AB2PRM 
MB 14  -  R3PRM+DB 1 -DA 1 08B*AB2PRM+DA1 0SA*AB3PRM 
WA11  -  -R1PRH+DA1+DA109B*AB3PRH+DA10SA*AB1PRM 
UB11  »  R3PRM+DB 1 -DA 1 OSB*AB 1PRM+DA10SA*AB3PRM 
WA12  -  R1PRH+DA1+DA108B*AB3PRM~DA108A*AB1PRM 
WB 12  -  R3P  RH+DB 1 4-DA  1 OSB  *AB  1 PRN+DA 1 08A*AB3PRM 
77  CONTINUE 

IF<TIME.  OT.  7.  000)00  TO  73 
RATE  -  -30.  /DEORAD 

CTHR1  -  -180.  /DEORAD  -  BTA11  +  BTA12 
CTHR2  -  -180.  /DEORAD  ♦  BTA22 
C 

C  GENERAL  EQUATIONS  FOR  INTEGRATION  DURING  ROTATION 
C 

TR2  -  TINE 
RD  ■  RATE  •  TR2 

D0HE1  -  -R1PRH  •  C08<RD)  -  R2PRM  •  SIN(RD) 

DOMES  -  -R1  *  COB(RD)  ♦  R3  •  BIN(RD) 

D08A1  •  A8F1«<AB1PRN*C08(RD>— AB2PRM*8IN(RD> ) 

♦AML 1 *AB3PRM* ( 1 .  +2.  *8IN(RD> )+£l*<AB3PRM*C08(RD) )+BAl 
D08B1  -  A8F1*AB3PRM-AM-1*< AB1PRM*C08 <RD)+AB2PRM*8IN< RD ) ) 
♦■B1-C1*(AB2PRM*C08(RD) )-El*A81PRM*< 1.  ♦SIN(RD) > 

D08A2  •  A8F2*(AB1*C08(RD)+AB3*BIN(RD> > +AML2*AB2+B A2 
♦E2*< AB2*SIN(RD) ) 

00888  -  A8F2*A82-AML8*<AB1«C08(RD)+AB3*SIN<RD) )+BB2 
.  -£2*<AB1*8IN<RD> >-C2*AB3*< 1.  -COB(RD) > 

00  TO  74 
73  CONTINUE 

TAa-08P3«aMEa+GHL2* ( CTHR2+R2*TR2 ) ♦DA2*TR2+DA208A*08A2 
♦DA208B*0B82 

TB2-08F2* ( CTHR2*R2*TR2 > -0ML2«OME2^DB2*TR2-DA2O8B*G8A2 
. ♦DA20SA*0S82 

TA1-08F1*0ME1+0ML1*<CTHR1+R3PRH*TR2)+DA1*TR2+DA108A*08A1 
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+DA10SB»08B1 

TB1-08F1*<CTHRH-R3PRM*TR2>-QML1*0NE1+DB1*TR2-DA10SB*0SA1 
+DA1 0SA*Q8B 1 

€♦♦♦♦♦♦♦»*♦»*♦♦*«♦»*♦»»*>****»»+****»■»♦»»»*♦♦*«**»**♦»♦♦*♦»♦♦■ 
THI8  SECTION  SIMULATES  THE  ROCKWELL-COLLINS 
PRE-FLIQHT  SENSOR  CALIBRATIONS.  IT  IS  ASSUMED 
THAT  ALL  BP IN-UP,  WAIT,  AND  ROTATIONAL  TIMES 
HAVE  BEEN  SATISFIED  AND  ACCELEROMETER  MEASUREMENTS 
ARE  AVAILABLE.  10  ITERATIONS  OP  THE  EQUATION 
SET  HILL  BE  U8ED  TO  MAKE  SURE  THAT  EXTREME 
ERR0R8  SETTLE  OUT  TO  OOOD  SOLUTIONS. 


THE  GRAVITY  AND  EARTH  RATE  COMPONENTS.  AS  MEASURED  BY 
AN  EXTERNALLY  LOCATED  I MU.  FOR  USE  IN  THE  CALIBRATION 
COMPUTATIONS ‘ARE  AS  FOLLOWS: 


AB1  -  QZEO 
AB2  -  OYEO 
AB3  -  OXEO 

AB1PRM  •  +DCMPPT *AB3  -  DCMPRL*AB2  +  AB1 
AB2PRM  •  -DCMPYA*AB3  +  AB2  ♦  DCMPRL*AB1 
AB3PRM  »  AB3  +  DCMPYA«AB2  -  DCMPPT#AB1 
R1  ■  OMEOZE 
R2  -  OMEOYE 
R3  *  OMEOXE 

R1PRM  -  +DCMPPT*R3  -  DCMPRL*R2  ♦  R1 
R2PRM  -  -DCMPYA*R3  ♦  R2  +DCMPRL*Rl 
R3PRM  -  R3  +  DCMPYA4R2  -  DCMPPT*Rl 
DO  110  I-l;10 

ASF2  -  <A21-A22-BTA22*AS3>/(2.  0+AB1 > 

AML2  *•  <0.  3*<B22-B21>+E2*AB3>/AB1 
BB2  -  .  5*  ( B2 1 +B  22 )  +E2*AB3— A8F2* AB2 

ASF1  -  (All-A12+<BTAll+BTA12)*AB2PRM)/<2.  0*(AB1PRM+E1*AB3PRM) > 
AML1  -  (0.  5*<B12-Bll)-E1*AB2PRM)/ABIPRM 
BB1  -  .  9*<BU+B12)+E1*AB1PRM-A8F1*AB3PRM 
AB2  -  <  B2 1 >AML2*AB 1 -BB2  > / ASF2 

AB3PRM  -  (B1 1+E1*AB2PRM+E1*AB1PRM+AML1«AB1PRM-BB1 )  /ASF1 

AB2PRM  -  -DCMP YA#AB3+AB2+DCMPRL*AB  1 

AB3  -  AB3PRM— DCMP YA*AB2PRM+DCMPP  T *AB 1 PRM 

BA2  -  .  3* ( A2 1 +A22+BT A22*AB3  > -AML2* AB2 

BA1  -  .  3*<AU«A12+<BTAll-BTA12>*AB2PRM)-AMLi*AB3PRM 

AB1  -  SORT < ABO* ABO  -  AB2*AB2  -  AB3*AB3 ) 

AB1PRM  •  SORT < ABO* ABO  -  AB2PRM*AB2PRM  -  AB3PRM*AB3PRM ) 

110  CONTINUE 

DA208A  •  (WA21-WA22+2.  *R1 > / (2.  *AB1 > 

DA20SB  »  ( WB22-WB2 1 ) / ( 2.  *AB1) 

DA2  -  0.  5*(WA21+WA22)-DA208B*AB2 
DB2  •  0.  S*<WB22+WB21)-R2-DA20SA*AB2 
DA108A  -  (UA 11 -HA 12+2.  *RlPRM)/<2.  *AB1PRM) 

DA108B  -  (WB12-WB1 1 ) / (2.  *AB1PRM) 

DAI  •  0.  9*(MA11+UA12>~DA108B*AB3PRM 

DB1  »  0.  S*<MBll+WB12)-R3PRM-Ai3PRM*DA10SA 

0ML2»<  <CTHR2+R2*TR2>*<TA2-DA2*TR2-DA208A*08A2-DA208»*«8B2> 


<-OME2*!T82-DB2*TR2-*DA208B*0SA2-DA20SA*«SB2) ) 

</<  <CTHR2+R2*TR2)**2*0ME2**2) 

0SF2- ( 0ME2*< TA2-DA2*TR2-DA2QSA*08A2-DA208B*Q8B2 ) ♦ ( CTHR2 
<+R2*TR2>* <  TB2-DB2*TR2*DA20SB*0SA2-DA20SA*0BB2 > ) 

</ < ( CTHR2+R2*TR2 >  **2+0ME2**2 ) 

©SF1— <  <TA1-DA1*TR2-DA1©8A*©9A1-DA10SB*0SB1 )*0ME1+(TB1 
<— DB 1 *TR2+DA 1 08B*GSA 1 —DA 1 08A*0S8 1 >*< CTHR1 +R3PRM*TR2) ) 

</  <  <  CTHR  1 *R  3PRM*TR2  >  **2+0ME  1  **2 ) 

0ML1-!  <TA1-DA1*TR2-DA1©SA*CSA1-DA10SB*0SB1 >*<CTHR1 
<+R3PRM*TR2)-0MEl*!TBl-DBl*TR2+DA10SB*08Al-DA108A*GSBl  > ) 
</!  !CTHR1+R3PRM*TR2>*»2+0ME1**2> 

PRINT*#  'CALIBRATED  TERMS  FROM  ROTATION' 

PRINT*#  * »» ******* »»♦»♦»♦♦♦»»♦»♦♦»»»»♦»*» * 

PRINT*,  'OML2.  OML1,  CSF1,  ©SF2-'#  ©ML2.  ©ML1. 08F1,  0SF2 
PRINT*#  '******************************** ' 

SDVX  »  AB3 
SDVV  -  ABB 
SDVZ  -  AB1 
C 

C ******  COMPUTE  CALIBRATED  ERROR  TERMS  FOR  ACCEL.  MODEL  EONS 
C 

C ******  WRITE  OUT  MCARLO  INPUT  VALUES 
WRITE (6. 137) 

WRITEI6. 134)  ASFX. ASFY. ASFZ 

WRITE (6, 139)  ABIABX. ABIASY,  ABIA8Z 

WRITE (6. 136)  EAA.EAAP 

ASFY  -  ASFY— A8F1+1.  O 

A3FZ  -  ASFZ— A8F2+1.  0 

ASFX  -  ASFY 

ABIA8X  -  ABIA8X-BB1 

ABIASY  -  ABIASY-(BAl+BB2)/2.  0 

ABIASZ  -  ABIA8Z—BAS 

EAAP  *  EAAP-AML1 

EAA  -  EAA-AML2 

WRITE (6,  138) 

WRITE!6. 134)  ASFX.  ASFY, ASFZ 

WRITE <6, 133)  ABIABX,  ABIASY,  ABIASZ 

WRITE (A. 136)  EAA.EAAP 

THT13  -  EAAP  ♦  EOAP 

THT31  -  EAA  *  EOA 

THT21  -  EAA 

BA11  -  1.  0 

BA12  -  P8II  *  TMT13 

BA13  -  -THTAI-THT12 

BA21  •  -PSII-THT23 

BA22  -10 

BA23  -  PHI I+THT21 

BA31  -  THTAH-THT32 

BA32  -  -PHI I-THT31 

BA33  -  1.  0 

134  FORMAT! 2X.  'ASFX.  ASFY.  ASFZ  -',3F12.  6) 

139  F0RMAT<2X.  'ABIASX,  ABIA8Y.  ABIASZ  -'#  3F12.  6> 

136  FORMAT (2X.  'EAA.EAAP  -',2F12.  6) 

137  FORMAT! 1H  ,  'MCARLO  INPUT  VALUES') 

138  FORMAT! 1H  #  'CALIBRATED  OUTPUT  ERRORS') 
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c««*« 

c«***«*  COMPUTE  CALIBRATED  ERROR  TERMS  FOR  OYRO  MODEL 

c###« 

08 FQ  -  OSFO  -  OS FI  >1.0 

OSFR  -  CSFR  -  0SF2  >1.0 

08FP  -  OSFO 

OB I ASP  •  OBJASP  -  DB1 

OBIASS  -  CBIA80  -  < DAl>DB2)/2.  0 

OBIASR  -  OBIASR  -  DA2 

AL1  -  AL1  -  DA208A 

AC1  -  AC1  -  DA20SB 

AL2  -  AL2  -  0A10SA 

AC  2  -  AC2  -  DA  10  SB 

EAW1  *  EAW1  -  0ML1 

EAW2  •  EAW2  -  0ML2 

PSI13  -  EAAP  >  EAW1 

P9IS1  -  EAA  >  EAUS 

B011  -  1.  0 

B012  •  P8II  >  PSI13 

B013  -  -THTAI  -  P8112 

B021  -  -P8II  -PSI23 

B022  -  1.0 

B023  -  PHI I  >  P8121 

B031  -  THTAI  >  P8I32 

B032  -  -PHI I  -  PBI31 

B033  -  1.  0 

STOP 

END 
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TABLE  A— 3.  EXAMPLE  RUNS 


[EXAMPLE  RUN  »  1 


INPUT  VALUES 


ABO 

AB1 

AB2 

32. 122547 

30.  664097 

-0.  057760  -9 

AB1PRM 

AB2PRM 

absprm 

30.  663145 

-0.  053737 

-9.  572191 

R1 

R2 

R3 

.  2331 14E-04  0.  353543E-04  0.  597B47E-04 

R1PRM 

R2PRM 

R3PRM 

.  222O9OE-04  0.  35346 IE -04  0.  5 9790 5E -04 

El 

E2 

0. 

174000E-02 

0.  174000E-02 

BTA11 

BTA12 

BTA22 

0. 

174000E— 02 

0.  174000E-02 

0.  174000E— 02 

0. 

DCMPRL 

DCMPPT 

DCMPYA 

0. 

100000E-03 

0.  100000E-03 

0.  loOOOOE— 03 

INPUT  VALUES 

TO  BE  CALIBRATED 

A8P2 

A8F1 

AML  2 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E— 03 

0. 

BA2 

BA1 

BB2 

0. 

488000E-01 

0.  488000E-01 

0. 48B000E-01 

0. 

AB1 

AB2 

AB3 

0. 

30664 1E+02 

-.  57760 5E -01 

-.  956912E+01 

AB1PRM 

AB2PRM 

ABSPRM 

0. 

30663 1E+02 

537372E-01 

-.  957219E+01 

DA20SA 

DA20SB 

DA108A 

0. 

145000E-02 

0.  145000E-02 

0.  145000E-02 

o. 

DA2 

DB2 

DAI 

0. 

970000E-05 

0.  97000QE-05 

0.  970000E-05 

o. 

CALIBRATION  ITERATION  OUTPUT 

A8F2 

A8F1 

AML2 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E-03 

0. 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E-03 

o. 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E-03 

0. 

0. 

lOOOOOE+Ol 

0. lOOOOOE+Ol 

0.  900000E-03 

o. 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0. 900000E-03 

0. 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E-03 

0. 

0. 

lOOOOOE+Ol 

0. lOOOOOE+Ol 

0.  900000E-03 

0. 

0. 

lOOOOOE+Ol 

0. lOOOOOE+Ol 

0.  900000E-03 

0. 

0. 

lOOOOOE+Ol 

0. lOOOOOE+Ol 

0. 900000E-03 

0 

0. 

lOOOOOE+Ol 

0.  lOOOOOE+Ol 

0.  900000E-03 

0. 

BA2 

BA1 

BB2 

0. 

488004E-01 

0.  487989E— 01 

0.  4BB000E-01 

0. 

0. 

488004E— 01 

0.  4B79S9E-01 

0.  4SB000E-01 

0. 

0. 

489004E-01 

0. 4879S9E— 01 

0. 4SB000E-01 

o. 

0. 

488004E-01 

0.  4S7989E-01 

0.  4SB000E-01 

o. 

0. 

488004 E -01 

0. 487989E— 01 

0. 48B000E-01 

0. 

0. 

4880O4E-O1 

0.  4B7989E-01 

0.  48SOOOE-01 

o. 

0. 

488004E-01 

0. 4S7989E— 01 

0.  4BS000E— 01 

0. 

0. 

4M004E-01 

0.  4S79B9E-01 

0. 4SS000E-01 

0. 

0. 

488004E-01 

0. 487989E-01 

0.  4SS0O0E-O1 

0. 

0. 

4SS004E— 01 

0.  4079B9E-O1 

0. 4SS000E-01 

0. 

AB1 

AB3 

AB3 

0. 

30664 IE +02 

57760 5E-01 

-. 956912E+01 

0. 

30664 IE +02 

577605E-01 

-.  956912E+01 

0.  30664 IE +03 

-.  57 760 5E -01 

-.  956912E+01 

0. 

30664 1E+02 

-.  577 60 5C -01 

-.  956912E+01 

0. 

306641E+02 

577 60 9E -01 

-.  956912E+01 

AB3 

569119 


BTA24 
1 74000E— 02 


AML1 

900000E-03 

BB1 

4B8000E-01 


DA10SB 
1 45000E-02 
DB1 

970000E-05 

AML1 

899990E— 03 
899990E-03 
899990E— 03 
899990E— 03 
899990E-03 
899990E— 03 
899990E-03 
899990E-03 
899990E-03 
899990E-03 
BB1 

4 8800 5E -01 
488014E— 01 
488014E-01 
488014E-01 
488014E-01 
488014E-O1 
488014E— 01 
488014E-01 
488014E— 01 
488014E-01 
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0.  306641 E+02  -. 377603E-01  -  9S6912E+01 
0.  30664 IE +02  -.  37760 3E -01  956912E+01 

0.  306641 E+02  377603E-01  9S6912E+01 

0.  30664 IE +02  3776036-01  936912E+01 

0.  30664 IE +02  377603E-01  936912E+01 

AB1PRM  AB2PRM  AB3PRM 

0.  30663 1E+02  -.  537372E-01  9372196+01 

0.  306631 E+02  337372E-01  957219E+01 

0.  30663 IE +02  337372E-01  937219E+01 

0.  30663 1E+02  337372E-01  -.  937219E+01 

0.  306631 E+02  337372E-01  93721 9E+01 

0.  306631 E+02  -.  S37372E-01  -.  937219E+01 

0.  306631 E+02  537372E-01  937219E+01 

0.  306631 E+02  537372E-01  937219E+01 

0.  306631 E+02  537372E-01  937219E+01 

0.  306631 E+02  537372E-01  937219E+01 

OUTPUT  FROM  OYRO  STATIONARY  MEASUREMENTS 

DA20SA  0A208B  DA1GSA  DA1CSB 

0.  14S000E-02  0.  1430006-02  0.  14S000E-02  O  143000E-02 
DA2  DB2  DAI  DB1 

0.  969954E-03  0.  970007E-05  0.  970139E-05  0  970066E-05 
0ML2  OMLl  CSF2  CSF1 

INPUT  VALUES  OP  OYRO  PARAMETERS 

0.  600000E-06  0.  600000E-06  0.  600000E-02  0.  600000E-02 
CALIBRATED  TERMS  FROM  ROTATION 
•*•****•*••«•+*•«••«•*•••••*••*••*•«*•****»•««•*•***• 
0ML2  OMLl  0SF2  GSF1 

O.  392792E— 06  O.  601671E-06  0.  600001E-02  O.  600000E-02 
*#*#*### ###*#*#*#*#***###+*##****•*****+*+ *********** 


EXAMPLE  RUN  •  2 


ABO  AB1  AB2  AB3 

32.122347  30.664097  -0.057760  -9.369119 

AB1PRM  AB2PRM  AB3PRM 

30.  661243  -0. 043691  -9.  376336 

R1  R2  R3 

.  2221 14E-04  0.  3S3343E-04  0.  597847E-04 
R1PRM  R2PRM  R3PRM 

.  22204 IE-04  O.  333297E-04  O.  398020E-04 
El  E2 

O.  174000E— 02  O.  174000E-02 

BTA11  BTA12  BTA22  BTA24 

0.  174000E-02  0.  174000E-02  0.  174000E-02  0.  174000E-02 
DCMPRL  DCMPPT  DCMPYA 

O.  300000E-03  0.  3000006-03  0.  300000E-03 
INPUT  VALUE8  TO  BE  CALIBRATED 

A8F2  A8F1  AML 2  AML1 

O.  1 OOOOOE+0 1  0.  100000E+01  0.  900000E-03  0.  900000E-03 
BA2  BA1  BB2  BB1 

0.  4880006-01  0. 488000E-01  0.  4880006-01  0.  488000E-01 
AB1  AB2  AB3 

0.  306641 E+02  37760 3E -01  9369126+01 

AB1PRM  AB2PRM  AB3PRM 

0.  306612E+02  -.  4369036-01  -.  9378346+01 

DA208A  DA208B  DA106A  DA10SB 

0.  1430006-02  0.  1430006-02  0.  1430006-02  0.  14S000E-02 
DA2  DB2  DAI  DBi 

0.  B70000E-03  0.  8700006-05  0. B70000E-03  0  8700006-03 
CALIBRATION  ITERATION  OUTPUT 
ASF2  A8P1 


AML2 


AML1 


0.  100000E+01  0.  100000E+01  0.  900000E-03  0.  900024E-03 
0.  lOOOOOE+Ol  0.  1000006+01  0.  900000E-03  0.  900029E-03 
0.  100000E+01  0.  lOOOOOE+Ol  0.  900000E-03  0.  900023E-03 
0.  lOOOOOE+Ol  O.  1000006+01  0.  900000E-03  0.  9OOO23E-03 
0.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  900000E-03  0.  900029E-03 
0.  1000006+01  0.  1000006+01  0.  900000E-03  0.  900023E-03 
O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  900000E-03  0.  900029E-03 
0.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  900000E-03  0  900029E-03 
0.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  900000E-03  0.  900029E-03 
O.  lOOOOOE+Ol  O.  1000006+01  0.  900000E-03  0.  900029E-03 
BAS  BAi  BBS  BB1 

O.  488OO4E-01  O.  487990E-01  0.  48B000E-01  0.  4B8009E-01 
0.  488004E-01  0.  487990E-01  0.  48B000E-01  0.  4B8024E-01 
0.  488004 E-Ol  0.  487990E-01  0.  48B000E-01  0.  488024E-01 
0.  4B8OO4E-01  0.  487990E-01  0.  488000E-01  0  488024E-01 
O.  488004 E -01  O.  487990E-01  0.  488000E-01  0.  488024E-01 
O.  488004E— 01  0.  4879906-01  0.  488000E-01  0  488024E-01 
O.  488004E-01  0.  487990E-01  0.  48B000E-01  0.  48B024E-01 
O.  488004E-01  0.  487990E-01  0.  48B000E-01  0  488024E-01 
0.  48800 4 E- 01  0.  4B799OE-01  0.  488000E-01  0  488024E-01 
0.  488004E— 01  0.  487990E-01  0.  48B000E-01  0.  488024E-01 
AB1  A82  AB3 

O.  306641E+02  977603E-01  936912E+01 

O.  30664 1E+02  977609E-01  996912E+01 

O.  30664 IE +02  9776096-01  -. 9969126+01 

O.  30664 IE +02  -.  9776096-01  996912E+01 

O.  30664 IE +02  -.  977609E-01  996912E+01 

O.  30664 IE +02  9776096-01  -.  9969126+01 

O.  30664 IE +02  9776096-01  936912E+01 

O.  30664 1E+02  -.  9776096-01  9969126+01 

O.  30664 1E+02  -.  9776096-01  996912E+01 

O.  3066416+02  9776096-01  996912E+01 

AB1PRM  AB2PRM  AB3PRH 

O.  306612E+02  4969096-01  937B34E+01 

O.  306612E+02  -.  4969096-01  -.  997834E+01 
O.  306612E+02  -.  4969096-01  997834E+01 

O.  306612E+02  -.  4 3690 9E -01  -.  9978346+01 
0.  306612E+02  -.  4969096-01  -.  997834E+01 
0.  306612E+02  -.  4969096-01  997834E+01 

O.  306612E+02  4969096-01  -.  9978346+01 

O.  3066126+02  496906E-01  -.  937834E+01 

0.  306612E+02  -.  496906E-01  -.9978346+01 
O.  3066126+02  -.  4969066-01  997B34E+01 

OUTPUT  FROM  OYRO  STATIONARY  MEASUREMENTS 

DA208A  DA208B  0A108A  DA10SB 

O.  1490006-02  0.  1490006-02  0.  1490006-02  O.  1490006-02 
0A2  DB2  DAI  DB1 

O.  870122E-03  0.  8701796-09  0.  B70321E-09  0.  870228E-09 
0ML2  0ML1  08F2  0SF1 

INPUT  VALUE8  OF  OYRO  PARAMETERS 

0.  4000006-06  0.  4000006-06  0.  4000006-02  0.  400000E-02 
CALIBRATED  TERMS  FROM  ROTATION 


0ML2  0ML1  OOFS  0SF1 

0.  39999 IE-06  0.  4048096-06  0.  4000016-02  0.  4000006-02 


ABO  AB1  AB2  AB3 

32.122947  30.664097  -0.097760  -9.369119 
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AB1PRM  AB2PRM  AB3PRM 

30.  661243  -0.  049691  -9.  576336 

R1  R2  R3 

.  2221 14E-04  0.  3535436-04  0.  997947E-04 
R1PRM  R2PRM  R3PRM 

.  22204 IE-04  0.  353297E-04  0.  398020E-04 
El  E2 

0.  174000E-02  0.  174000E-02 

BTA11  BTA12  BTA22  BTA24 

0.  174000E-02  0.  174000E-02  0.  174000E-02  O  174000E-02 
DCMPRL  DCMPPT  DCMPYA 

0.  300000E-03  0.  300000E-03  0.  300000E-03 
INPUT  VALUES  TO  BE  CALIBRATED 

ASF 2  A8F1  AML 2  AML1 

O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  400000E-04  0  400000E-04 
BA2  BA1  BB2  BB1 

O.  322000E-01  0.  322000E-01  0.  322000E-01  0  322000E-01 
AB1  AB2  AB3 

O.  306641E+02  577605E-01  936912E+01 

AB1PRM  AB2PRM  AB3PRM 

0.  306612E+-.’  4 5690 5E -01  -  957B34E+01 

DA20SA  DA208B  DA10SA  DA1CSB 

O.  145000E-02  0.  149000E-02  0.  14S000E-02  0.  145000E-02 
DA2  DB2  DAI  DB1 

0.  970000E-05  0.  970000E-05  0.  970000E-05  O  9700O0E-0S 
CALIBRATION  ITERATION  OUTPUT 

ASF2  ASF1  AML2  AML1 

O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  O.  400000E-04  0  400105E-04 
O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  399997E-04  0  400109E-04 
0.  lOOOOOE+Ol  O.  lOOOOOE+Ol  0.  399997E-04  O.  400105E-04 
O.  lOOOOOE+Ol  O.  lOOOOOE+Ol  O.  399997E-04  O  400105E-04 
O.  lOOOOOE+Ol  O.  lOOOOOE+Ol  0.  399997E-04  O.  400109E-04 
0.  lOOOOOE+Ol  O.  lOOOOOE+Ol  0.  399997E-04  O.  400105E-04 
O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  399997E-04  O  400105E-04 
O.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  399997E-04  O.  400105E-04 
0.  lOOOOOE+Ol  0.  lOOOOOE+Ol  0.  399997E-04  O.  400105E-04 
0.  lOOOOOE+Ol  O.  lOOOOOE+Ol  0.  399997E-04  O.  400105E-04 
BA2  BA1  BB2  BB1 

0.  322007E-01  0.  321V97E  01  0.  322000E-01  0  321999E-01 

0.  322007E— 01  0.  321997E-01  0.  322000E-01  0  32201 BE -01 

O.  322007E— 01  O.  321 99 7E -01  0.  322000E-01  O.  32201 BE -01 

O.  322007 E -01  0.  321997E-01  0.  322000E-01  O.  322016E-01 

O.  322007E-01  0.  321997E-01  0.  322000E-01  0  322016E-01 

O.  322007E-01  0.  32 1 99 7 E -01  0. 322000E-01  O.  32201 BE -01 

0.  322007 E -01  0.  321997E-01  0.  322000E-01  0.  32201 BE -01 

0.  322007E-01  0.  321997E-01  0.  322000E-01  0  32201 BE -01 

O.  322007 E -01  0.  321997E-01  0.  322000E-01  0  32201BE-01 

0.  322007 E -01  0. 321997E-01  0.  322000E-01  0.  32201 BE -01 

AB1  AB2  AB3 

0.  30664 1E+02  -.  37760 3E -01  99691 2E+01 

0.  30664 IE +02  577605E-01  9969126+01 

0.  30664 1E+02  577609E-01  996912E+01 

0.  30664 IE +02  977609E-01  996912E+01 

0.  30664 1 E+02  37760 SE -01  996912E+01 

0.  306641 E+02  -.  577609E-01  996912E+01 

0.  306641 E+02  977605E-01  996912E+01 

0.  30664 IE +02  977609E-01  996912E+01 

0.  306641 E+02  377 60 3E -01  996912E+01 

0.  306641 E+02  377603E-01  996912E+01 

AB1PRM  AB2PRM  AB3PRM 

0.  306612E+02  496909E-01  997834E+01 
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0.  306612E+02  45690SE-01  957B34E+01 

0.  3O6612E+02  4 5690 5E -01  957B34E+01 

0.  306612E+02  456905E-01  957B34E+01 

0.  306612E+02  456905E-01  957834E+01 

0.  306612E+02  45690 5E-01  957B34E+01 

0.  306612E+02  45690 5E -01  957B34E+01 

0.  3O6612E+02  456906E-01  957834E+01 

0.  3O6612E+02  -  456906E-01  957834E+01 

0.  306612E+02  4 5690 6E -01  -  957834E+01 

OUTPUT  FROM  OYRO  STATIONARY  MEASUREMENTS 

DA20SA  DA20SB  DA10SA  DA1SSB 

0.  145000E-02  0.  145000E-02  0.  145000E-02  0.  145000E-02 
DA2  DB2  DAI  DB1 

0.  969959E-05  O.  97001 3E-05  0.  970252E-05  0.  970438E-05 
QM  9  OML1  0RF9  occi 

INPUT  VALUES  OF  OYRO  PARAMETERS 

0.  500000E-06  0.  500000E-06  0.  900000E-02  0.  500000E-02 
CALIBRATED  TERM8  FROM  ROTATION 
•***•*•**«*••«*«***'»•******••«**•*••******«•«••**•*' 
0ML2  0ML1  0SF2  CSF1 

O.  48874 5E-06  0.  50052 IE-06  0.  50000 IE-02  0.  500000E-02 


DISTRIBUTION 


No.  of 
Copies 

DR SMI -R  1 

DRSMI-LP ,  Mr.  Voigt  1 

PRSMI-RG,  Dr.  Yates  1 

DRSMI-RGN,  Mr.  McLean  6 

DRSMI-RGT,  Mr.  Jordan  1 

DRSMI-RPT  (Record  Copy)  1 

DRSMI-RPR  15 

US  Army  Material  System  Analysis  Activity  1 

ATTN:  DRXSY-MP 

Aberdeen  Proving  Ground,  MD  21005 
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